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METHOD AND DEVICES FOR LASER tissue. Furthermore, backscattering does not account for the 

INDUCED FLUORESCENCE ATTENUATION effects of absorption and scattering suffered by the intrinsic 

SPECTROSCOPY fluorescence prior to its measurement. 

r> A^^D/^TTivTn ac titc TM^NmnM In biological LIF techniques, lasers are used to cause 

BACKGROUND OF THE INVENTION 5 fluorophores fa the sample t0 emit fluorescence . The main 

The present invention is directed to methods and devices fluorophores in normal biological tissue are tryptophan, 

for determining a spectroscopic characteristic of a sample collagen and elastin. Other fluorophores, such as NAD 

utilizing laser induced fluorescence attenuation spectros- (nicotinamide adenine dinucleotide ) and FAD (flavin 

copy ("LIFAS" ). More particularly, the invention is directed adenine dinucleotide), are also normally present, but at 

to methods and devices for measurement of the wavelength- 10 much lower concentrations. Under certain conditions, the 

dependent attenuation of the sample and subsequent resto- contribution of the various fluorophores to the LIF can 

ration of the intrinsic laser induced fluorescence ("LIF') for change. For example, during an ischemic or hypoxic event, 

physiological monitoring, biological tissue characterization the tissue is deprived of oxygen and anaerobic respiration 

and biochemical analysis. takes place. Consequently, the weak fluorophore NAD will 

Conventionally, samples have been characterized by 15 be converted into the strong fluorophore reduced nicotina- 

determining the attenuation and laser.induced fluorescence mide ademne dinucleotide CNADfT). As a result, the LIF 

("LIF'). Once the attenuation and LIF of a sample have been collected from the sample will reflect an increased contn- 

determined, these spectroscopic properties can be utilized to butlon from NAE> H. This change is typically observed as a 

determine a physical or physiological property of the 0A rise in the intensity of the LIF spectrum in the region of peak 

sample. For example, the attenuation of a sample can be "° NADH emission, from about 470 to 490 nm. Thus, the 

used to determine the concentrations of mixture components metabolic state of the tissue can be determined by measuring 

orturbidity of afluid. Similarly, the LIF of a sample has been the relative chan S e in LIF intensity at the wavelength of 

used in fields such as analytical chemistry, environmental P eak NADH emis s ion as compared to the relative change of 

monitoring, industrial inspection and medical diagnosis. In 0 the LIF ^tensity of fluorophores that are normally present in 

the medical field, for example, LIF spectroscopic techniques me tissue > such ^ eIastin or collagen, 

have been used for tissue characterization, malignant tumor Unfortunately, the LIF of biological tissue is heavily 

identification, atherosclerotic plaque diagnosis, metabolism modulated in the 390-450 nm range by the peak absorption 

evaluation, and the like. of the main tissue chromophore, hemoglobin. Thus, 

Conventionally, the attenuation or the absorption of a 30 although the intrinsic LIF spectra of normal tissue should 
sample is determined by placing the sample between a light approximately resemble that of the pure fluorophore corn- 
source and a detector and measuring any reduction in the ponents of the tissue, the measured LIF spectra has a valley 
intensity of the light as it passes through the sample. In order in * e 4 00-450 nm region that is associated with the hemo- 
to obtain a measurement having an acceptable signal-to- S lobin absorption. As a result, the measured LIF spectrum 
noise ratio using these conventional techniques, it is impor- ^ tissue appears to have a double peak instead of the single 
tant to transmit the incident light with sufficient intensity. " P eak spectrum associated with the pure fluorophores of the 
Thus, the thickness of the sample and the wavelength of the tissue * Hence, the LIF spectrum of normal tissue begins to 
incident light are important factors affecting the reliability of resemble the spectrum associated with tissue suffering 
the resulting measurement. Moreover, because it is neces- hypoxia or ischemia, which makes it more difficult to 
sary to place the sample between the light source and the 40 identif Y tissue abnormalities. 

detector, it is difficult to perform attenuation measurements Furthermore, since the optical properties of biological 

on certain types of samples, such as living tissue. tissue are influenced by its hemoglobin concentration, the 

More recently, fiber optic techniques have been developed measured LIF will vary with the level of blood perfusion 

for the measurement of attenuation in which an optical fiber throughout the cardiac cycle. The LIF of contractile tissue, 

is used to guide the incident light to illuminate a sample 45 sucn 35 & e myocardium, is highly dependent on its state of 

inside a small chamber at the tip of a probe. A reflector is contraction. Contraction increases the concentration of the 

placed on the opposite side of the chamber to reflect the fluorophore NADH and, hence, its contribution to tissue 

incident light into a second fiber which is associated with a fluorescence. During contraction, however, blood is pumped 

detector. Unfortunately, these techniques find limited appli- out of ^e tissue, thereby reducing its hemoglobin concen- 

cation with materials, such as fluids, that can readily pass 50 Nation and, hence, light attenuation. Thus, a contracted 

into the chamber in the tip of the probe. myocardium retains less blood (i.e., a lower hemoglobin 

Conventionally, laser induced fluorescence CLIFS") concentration) and, therefore, exhibits lower light 

spectroscopic techniques utilize various optical configura- absorption, than when relaxed. 

tions in which a laser is directed at a sample using an optical Organs such as the brain, heart and kidney are the most 

fiber and the LIF from the sample is collected using a second 55 sensitive to oxygen deficiency and can suffer permanent 

optical fiber. Alternatively, the same fiber can be used for damage following an ischemic or hypoxic event. During 

excitation of the sample and the collection of LIF. In either open-heart surgery, for example, continuous monitoring of 

case, the LIF collected by the fiber is modulated by the kidney perfusion, i.e., ischemia, is required. Similarly, since 

sample, e.g. by the wavelength-dependent absorption and me success of a transplantation surgery is highly dependent 

scattering of constituents of the sample. Therefore, existing 60 on me level of organ perfusion at harvest and during 

LIFS methods are limited by the fact that the "intrinsic" or preservation of the tissue, continuous monitoring of the 

"true" fluorescence of the sample's fluorophores cannot be or gan is typically required. 

determined. Recent reports have suggested measuring the Ischemia and hypoxia are both conditions that deprive 

diffuse reflectance spectrum of the tissue as a means for tissue of oxygen, leading to anaerobic metabolism and the 

correcting the intrinsic LIF using Monte-Carlo mathematical 65 accumulation of the metabolic coenzyme NADH. The coen- 

formulations. However, such correction methods are criti- zyme NADH is a fluorescent molecule. Therefore, ischemia 

cally dependent on the backscanering characteristics of the and hypoxia can be indirectly detected using LIF techniques 
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by sensing increased concentrations of NADH and inter- 
preting its elevation as a sign of oxygen deficiency. A 
common indicator of oxygen deficiency is the ratio between 
the LIF intensity at wavelengths associated with the peak 
fluorescence emission of NADH, collagen and elastin. 5 
However, such methods have not been practically applied 
for the detection of ischemia because of several complica- 
tions. First, these methods cannot determine whether the 
elevated NADH concentration is caused by ischemia, 
hypoxia or hypermetabolism. Second, scarred or fibrosed 10 
tissue would be detected as normal because of their low 
NADH concentration. Finally, the indicator ratios are cal- 
culated by normalizing the intensity of NADH peak fluo- 
rescence by that of collagen or elastin. Although the fluo- 
rescence of the structural proteins elastin and collagen does 15 
not vary with tissue oxygenation, their fluorescence vary 
with the site of measurement. 

The intrinsic LIF of biological tissue has a single peak 
spectral profile resembling that of the pure structural pro- 
teins which are primarily elastin and collagen. However, this 
single peak profile is modulated by the attenuation of 
hemoglobin, especially about 410 nm, to yield a distorted, 
double-peak spectrum. The second spurious peak at about 
470 nm results from the spectral valley created by heavy 
hemoglobin absorption at about 410 nm. Unfortunately, this 
second peak overlaps and obscures the peak fluorescence 
emission of NADH at about 470-490 nm, which may 
impede ischemia detection techniques that are based on 
NADH concentration. 

The methods and devices of the present invention will 30 
enable measurement of several LIF parameters that directly 
relate to the state of tissue perfusion and/or hemoglobin 
oxygenation. This offers superior accuracy over the cur- 
rently available LIF techniques which detect ischemia and ^ 
hypoxia indirectly by sensing an upsurge in NADH fluo- 
rescence. In particular, the methods and devices of the 
present invention can directly indicate the presence and the 
level of hypoxia or ischemia. Furthermore, the methods of 
the present invention utilize novel criteria for discrimination ^ 
between normal and ischemic or hypoxic tissue. 

It is readily apparent that existing LIFS techniques in 
medicine and industry are affected by modulations of the 
mtrinsic LIF. Thus, there is need for methods and devices for 
determining the intrinsic laser induced fluorescence while 45 
eliminating or minimizing modulating influences in order to 
provide more accurate and complete characterization of a 
sample. In particular, there is a need for method and devices 
for characterization of biological tissue. The present inven- 
tion addresses these and other needs. 5( 

SUMMARY OF THE INVENTION 

The present invention is embodied in a system and related 
method for laser-induced fluorescence attenuation spectros- 
copy ("LIFAS") in which the attenuation and intrinsic LIF of 55 
a sample can be determined. The LIFAS system employs a 
source, a first sensor, a second sensor and a processor. The 
source, preferably a laser, emits light to irradiate a sample 
volume in a sample so that the sample volume produces 
return light, which preferably includes laser induced fluo- 60 
rescence. The first sensor monitors the return light at a first 
distance from the sample volume and generates a plurality of 
signals representing the intensity of the return light in 
predetermined wavelength bands. The second sensor moni- 
tors the return light at a second distance from the sample 65 
volume and generates a plurality of signals representing the 
intensity of the return light, preferably over the same wave- 
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length bands. Where the first and second distances are 
different, the processor can be used to determine the 
wavelength-dependent attenuation of the sample using the 
signals of both sensors. The measured attenuation will 
typically reflect the effects of both absorption and scattering 
by the sample. For maximum signal-to-noise ratio, it is 
preferred that the first and second detectors monitor return 
light from the sample volume at a location in proximity to 
the sample volume. 

Once the wavelength-dependent light attenuation is 
determined, the attenuation can be used to restore the 
intrinsic LIF of the fluorophores of the sample. Once the 
attenuation and intrinsic fluorescence have been determined, 
these factors can be used in a variety of applications, 
including medical diagnosis. In the preferred embodiment, 
the attenuation and fluorescence are used for the detection of 
ischemia and hypoxia of biological tissue. Ischemia and 
hypoxia are detected by monitoring spectral changes in the 
intrinsic LIF to detect a rise in the concentration of NADH 
20 in the tissue. Since ischemia is caused by a decrease in the 
blood content of the tissue, ischemic tissue will have lower 
concentrations of hemoglobin. Consequently, in ischemic 
tissue, absorption by hemoglobin, especially in the 390-450 
nm spectral band, will be reduced. Changes in the wave- 
25 length of peak light transmission, especially in the 450-500 
nm band, can also be used to detect hemoglobin oxygenation 
and, hence, the presence of hypoxia or ischemia. Finally, 
multicriteria associative memories (MAM) can be used as an 
attenuation or fluorescence spectral classifier. 

Other features and advantages of the present invention 
will be apparent from the following description of the 
preferred embodiments, taken in conjunction with the 
accompanying drawings, which illustrate, by way of 
example, the principles of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a schematic diagram of a spectroscopic system 
in accordance with the invention; 

FIG. 2(a) is a schematic diagram of a Laser Induced 
Fluorescence Attenuation Spectroscopy system O'LIFAS") 
in accordance with the present invention having an excita- 
tion waveguide and a first and second collection 
waveguides; 

FIG. 2(b) is a schematic diagram of an alternative con- 
figuration of the collection waveguides of FIG. 2(a); 

FIG. 2(c) is a schematic diagram of yet another alternative 
configuration of the collection waveguides of FIG. 2(a); 

FIG. 3(a) is a schematic diagram of a LIFAS system in 
accordance with the present invention having a first 
excitation-collection waveguide and a second collection 
only waveguide; 

FIG. 3(&) is a schematic diagram of an alternative con- 
figuration of the collection waveguides of FIG. Ma); 

FIG. 4 is schematic diagram of a LIFAS system for 
biomedical applications in which excitation radiation is used 
to induce fluorescence of tissue; 

FIGS. 5(a) and (b) are end perspective views of a LIFAS 
probe having an optical fiber and an optical fiber bundle as 
the central waveguide; 

FIGS. 6(a) and (b) are partial perspective views of an 
electrical subsystem used to synchronize a LIFAS system 
with the electrical activity of tissue using electrodes inte- 
grated into the optical probe; 

FIGS. 1(a) and (b) are schematic diagrams showing the 
relative positions of the action potential and electrodes 
integrated into the optical probe; and 



6,12 

5 

FIGS. 7(c) and (d) are the corresponding action potential 
timing diagrams; 

FIG. 8 is a perspective view of a probe incorporating a 
visible light source for probe positioning and a circular LED 
array for indicating the direction of contraction propagation; 

FIGS. 9(flHc) are partial perspective views of three 
LIFAS probes having waveguides arranged in different 
geometrical configurations adapted for different applica- 
tions; 

FIG. 10 is a schematic diagram of a biomedical LIFAS 
system in accordance with the invention in which a single 
optical detector is used to measure the return light from two 
collection fibers; 

FIG. 11 is a block diagram of the method for determining 
the attenuation and the intrinsic fluorescence of a sample; 

FIG. 12 is a graph of the attenuation spectrum measured 
from normal and ischemic renal cortex of rabbit kidney 
using LIFAS devices and methods in accordance with the 
invention; 

FIG. 13 is a scatter plot of the fluorescence intensity at 
480 nm acquired from hypoxic (+), normal (o) and hyper- 
oxic (x) tissue through excitation-collection and collection- 
only waveguides; 

FIG. 14 is a graph of the intensity l co (X) c of the signal 
measured by a collection-only waveguide from normal, 
hyperoxic and hypoxic tissue; 

FIG. 15(a) is a graph of the modulated LIF spectra of 
normal and ischemic kidney, and 

FIG. 15(b) is a graph of the intrinsic LIF spectra of normal 
and ischemic kidney; 

FIGS. 16(a) and (b) are graphs of typical modulated LIF 
spectra of normal and ischemic kidney, respectively, depict- 
ing the difference in the symmetry of the main lobe; and 

FIGS. Yl(a)Ab) are graphs of the mean LIF intensity of 
live kidney and heart tissue and 

FIGS. \l(c)Ad) are the corresponding LIFAS-derived 
attenuation spectra. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

As shown in the exemplary drawings and, in particular, 
FIG. 1 thereof, the present invention is embodied in a 
spectroscopic system 110 and related method for measuring 
the attenuation and/or optical rotation caused by a sample 
114. The spectroscopic system 110 shown in FIG. 1, 
includes a source 111 that produces radiation 112, which is 
directed at a sample volume 113 within the sample 114. The 
source 111 is preferably a laser and the radiation 112 is 
preferably monochromatic ultraviolet (UV), visible or infra- 
red (IR) radiation 112. Nevertheless, the source 111 can be 
selected to produce other types of radiation, such as broad- 
band or polarized radiation, useful for a particular applica- 
tion. The spectroscopic system 110 further includes a first 
sensor 116 displaced by a distance Zl from the sample 
volume 113, and a second sensor 118 displaced by a distance 
72 from the sample volume 113, which are used to monitor 
return radiation 120 from the sample 114. In the preferred 
embodiment, the return radiation 120 will include fluores- 
cence of the fluorophores of the sample 114. 

The sensors 116 and 118 each include an aperture 116a 
and 118a, respectively, through which the return radiation 
120 is observed and collected. The sensors 116 and 118 
typically can only observe the events occurring within the 
volume defined by a solid angle emerging from the aperture 
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of the sensor. This volume is known as the numerical 
aperture or the field of view of the sensor. Thus, first and 
second sensors 116 and 118 include numerical apertures 117 
and 119, respectively, that are adapted to receive a portion of 
5 the return radiation 120 from the sample 114. The sensors 
116 and 118 can be tilted at different angles O, and 8 2 and 
selected to have different numerical apertures 117 and 119 so 
that the sensors 116 and 118 can selectively monitor return 
radiation 120 from all or part of the sample volume 113. 
10 The sensors 116 and 118 each generate a signal 124a and 
125a, respectively, representing the intensity I rl (X) and 
l c ^(X) and/or polarization of the return radiation 120. 
Preferably, the signals 124a and 125a are generated at a 
plurality of wavelengths within predetermined wavelength 
15 bands. The signal 124a from the first sensor 116 and the 
signal 125a from the second sensor 118 are communicated 
by signal paths 124 and 125, respectively, to a processor 
123. As discussed in additional detail, below, the processor 
123 processes the signals 124a and 125a to determine a 
20 spectral characteristic of interest, such as the attenuation 
and/or optical rotation of the return radiation 120 caused by 
the sample 114. 

To determine the attenuation of the return radiation 120 
caused by the sample, the sensors 116 and 118 are positioned 
25 at non-equivalent distances Zl and Z2 from the sample 
volume 113. Preferably, the first and second sensors 116 and 
118 will be positioned so that one of the sensors will be in 
the immediate proximity of the sample volume 113 while the 
other sensor is positioned at a location adjacent to the first 
30 sensor, but displaced slightly further from the sample vol- 
ume 113. It has been found that this arrangement advanta- 
geously improves the signal-to-noise ratio and, hence, mea- 
surement accuracy of the spectroscopic system 110. 
The attenuation of the return radiation 120 is dependent 
35 on the distance traveled by the return radiation 120 through 
the sample 114. Hence, the signal 124a measured by the first 
sensor 116 and the signal 125a measured by the second 
sensor 118 will suffer different attenuation where the dis- 
tances Zl and Z2 from the sample volume 113 are not 
40 equivalent. The processor 123 can be used to process the 
signals 124a and 125a to determine the wavelength - 
dependent attenuation of the sample 114. The attenuation 
determined by the spectroscopic system 110 will reflect the 
effects of light absorption and scattering by the sample 114. 
4 Where the radiation 112 causes fluorophores in the sample 
114 to fluorescence, the attenuation of the sample 114 can be 
used to determine the intrinsic fluorescence 129 of the 
fluorophores in the sample volume 113. The methods by 
5Q which the attenuation and intrinsic LIF are determined are 
discussed in detail below. 

The source 111, shown in FIG. 1, can be any radiation 
source, such as a laser or a lamp emitting excitation radiation 
112 at a wavelength capable of interacting with the sample 
55 114. The sensors 116 and 118 will preferably include a 
detector, which can be as simple as individual light-sensitive 
diodes with appropriate band pass filters, or more compli- 
cated sensors such as an optical spectrum analyzer. 
Preferably, the sensor is a suitable spectrograph or spec- 
50 trometer equipped with a suitable sensor. Alternatively, the 
sensor can include a multispectral CCD camera. 

In the preferred embodiment, the source. Ill is an XeCl 
excimer laser emitting monochromatic ultraviolet radiation 
at 308 nm. Each of the sensors 116 and 118 is a spectrograph 
65 (Model FF250, ARIES Inc., Concorde, Mass.) associated 
with a 1024 element intensified photo diode array (PDA) to 
image the resolved spectrum. Each PDA is connected to an 
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optical multichannel analyzer (OMA III, EG&G Princeton 
Applied Research Corporation, Princeton, N.J.). The OMA 
reads the light spectrum imaged by the PDA to produce the 
signals 124a and 125a. A 335 nm longpass filter (Schott 
WG335) is placed in front of the aperture of each of the 
sensors 116 and 118 to remove any backscattered 308 nm 
excitation radiation. The processor 123 is a personal com- 
puter that is networked to each of the sensors 116 and 118 
via the signal paths 124 and 125 to receive the signals 124a 
and 125a, respectively. 

In other applications, such as measuring the optical rota- 
tion or anisotropy of a sample, it will be desirable to position 
the first and second sensors 116 and 118 at equivalent 
distances Zl and 21 from the sample volume 113. Optical 
rotation is observed in optically active materials. An opti- 
cally active material is characterized by a lack of symmetry 
in its molecular or crystalline structure which causes a 
rotation of the plane of polarization in incident-plane- 
polarized radiation. The extent to which the plane of polar- 
ization is rotated will typically vary from one optically 
active material to another. In addition, the extent of rotation 
may depend on the number of molecules in the path of the 
radiation, the wavelength of the radiation, and the tempera- 
ture of the material. Anisotropy, on the other hand, is caused 
by a selective response of the sample to radiation travelling 
in different directions. In certain biological samples, for 
example, anisotropy can result from microscopic inhomo- 
geneity in the tissue structure. Use of equivalent distances 
Zl and Z2 will minimize the contribution of the path-length- 
dependent attenuation to the overall attenuation measured by 
the system 110. Therefore, the optical rotation or anisotropy 
of the sample can be determined more easily. 

The LIFA system 110 shown in FIG. 1 can be adapted to 
measure optical rotation by utilizing first and second sensors 
116 and 118 which are adapted to measure the polarization- 
dependent intensity of the return radiation. In particular, the 
first sensor 116 can be adapted to measure a first optical 
angle <j>, (not shown) yielding the maximum intensity of the 
first portion of the return radiation 120. Meanwhile, the 
second sensor 118 can be adapted to measure a second 
optical angle <j) 2 (not shown) yielding the maximum intensity 
of the second portion of the return radiation 120. The optical 
rotation of the return light 120 will be given by the differ- 
ence between the angles 0, and provided that the system 
has been properly calibrated. The same sensor is used for 
determination of the attenuation of the sample will find 
utility for monitoring the optical rotation. These sensors can 
be adapted to such measurements by placing a rotatable 
polarizing filter at each of the apertures 116a and 118a of the 
first and second sensors 116 and 118, respectively. 

An alternative embodiment of the present invention is 
shown in FIG. 2(a). In this embodiment, the spectroscopic 
system 110 of FIG. 1 has been adapted as a Laser Induced 
Fluorescence Attenuation Spectroscopy ("LIFAS") system 
210. In the LIFAS system 210, a source 211 emits laser 
radiation 212 at an intensity and a wavelength capable of 
inducing fluorescence of the sample 214. The laser radiation 

212 is transmitted through an excitation waveguide 215 to a 
sample volume 213 within the sample 214 where the laser 
radiation 212 excites local fluorophores in the sample vol- 
ume 213 to emit intrinsic fluorescence 229. The intrinsic 
fluorescence 229 of the fluorophores in the sample volume 

213 is modulated, for example, by the absorption and 
scattering of the local chromophores (not shown) and scat- 
tered (not shown) present in the sample 214, respectively. 

Two collection waveguides 221 and 222 having apertures 
221a and 222a and numerical apertures 217 and 219 are 
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disposed about the excitation waveguide 215. As shown in 
FIG. 2(b), the collection waveguides 221 and 222 may be 
positioned such mat their apertures 221a and 222b are 
laterally displaced from the excitation waveguide 215 by the 
distances xT and x2' and/or axially . displaced from the 
aperture 215a of the excitation waveguide 215 by the 
distances y 1 ' and y2\ Preferably, the lateral distances xl' and 
x2' are small or zero. Meanwhile, the axial distances y 1' and 
y2' are preferably selected to be small and unequal, while 
one of the axial distances is preferably zero. Also, it is 
desirable to position the apertures 215a, 221a and 222a of 
the waveguides 215, 221 and 222 in close proximity to or in 
contact with the sample 214 during use of the LIFAS system 
210, as shown in FIG. 2(a). 

An alternative configuration of the collection waveguides 
is shown in FIG. 2 (c). In this configuration, the axial 
distances yV and y2* are both zero. Meanwhile, the lateral 
distances xl' and x2' are selected to be small and unequal. 
The apertures 240a and 242a of the waveguides 240 and 242 
are oblique such that the fields of view 244 and 246 are 
directed toward the sample volume 213 and to advanta- 
geously encompass the bulk of the sample volume 213. As 
shown in FIG. 2(c), the distances zl' and z2' are directed 
along the respective axes of the numerical apertures 244 and 
246 and extend from the surface of the apertures 240a and 
242a to the point of intersection of the axes of the numerical 
apertures 244 and 246, respectively. 

The collection waveguide 221 collects a first portion 228 
of the return light 220 and transmits the first portion 228 of 
the return light 220 to a first sensor 216. The first sensor 216 
generates a first signal 224a. representing the intensity I c ,(^) 
of the first portion 228 of the return light 220 at a plurality 
of wavelengths within predetermined wavelength bands. 
Similarly, the collection waveguide 222 collects a second 
portion 230 of the return light 220 and transmits the second 
portion 230 of the return light 220 to the second sensor 218. 
The second sensor 218 generates a second signal 225a 
representing the intensity I c2 (X) of the second portion 230 of 
the return light 220 at a plurality of wavelengths preferably 
within the same wavelengths bands monitored by the first 
sensor 216. 

The processor 223 then processes the first and second 
signals 224a and 225a to determine the wave length - 
dependent attenuation of the sample 214, As discussed in 
detail, below, once the attenuation is known, either of the 
signals 224a or 225a can be processed to minimize the 
effects of attenuation and determine the intrinsic LIF 229 of 
the fluorophores in the sample volume 213. Alternatively, if 
the sensors 216 and 218 are selected to monitor the polar- 
ization of the first portion 228 and the second portion 230 of 
the return light 220, the processor 223 can be used to 
determine the optical rotation of the return light 220 caused 
by the sample 214. 

Preferably, the intensity and wavelength of the laser 
radiation 212 should be sufficient to create a sample volume 
213 that is large enough to have some overlap with the 
numerical apertures 217 and 219 of the collection 
waveguides 221 and 222, respectively. The numerical aper- 
tures 217 and 219 can be selected to encompass all or parts 
of the sample volume 213. Furthermore, longpass filters are 
placed before the inputs of the sensors 216 and 218 to 
selectively block backscattered excitation radiation. 

In the preferred embodiment, the source 211 is an Xe — CI 
excimer laser emitting 308 nm ultraviolet excitation radia- 
tion 212. The waveguides 215. 221 and 222 are advanta- 
geously comprised of optical fibers or optical fiber bundles. 
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which can be integrated into a probe for ease of use and 
durability. The optical bundle and fibers are made of fused 
silica which is transparent to ultraviolet radiation. 
Specifically, a 1.4 mm diameter optical bundle is used as the 
excitation waveguide 215 while the collection waveguides 5 
221 and 222 are 0.4 mm optical fibers disposed about the 
excitation bundle. In the preferred probe embodiment, the 
lateral range xl* and x2* are equal to zero while the axial 
distance yl' is about 0.6 mm and the axial distance yT is 
about zero. Furthermore, the probe may be configured J0 
within, for example, the shaft of a hypodermic needle, so 
that the in-vivo physiological or pathological properties of 
biological tissue can be determined. 

In the preferred embodiment, each of the sensors 216 and 
218 is a spectrograph (Model FF250, ARIES Inc., Concorde, i5 
Mass.) associated with a 1024 element intensified photo 
diode array (PDA). Each PDA is connected to an optical 
multichannel analyzer (OMA III, EG&G Princeton Applied 
Research Corporation, Princeton, N.J.) which measures the 
intensity of the light spectrum imaged by the PDA and 20 
produces the signals 224a and 225a. The processor 223 is a 
personal computer that is networked to each of the sensors 
216 and 218 via the signal paths 224 and 225 to receive the 
signals 22Aa and 225a, respectively. 

The method of determining the attenuation and intrinsic 2 5 
fluorescence of a sample is better understood with reference 
to the block diagram shown in FIG. 11, which illustrates the 
main steps involved in determining the wavelength- 
dependent attenuation coefficient a(X) and the intrinsic 
fluorescence I^X) of a sample in accordance with the current 30 
invention. In the embodiment shown in FIG. 2(a), the first 
and second portions 228 and 230 of the return light 220 
experience different attenuation effects due to the unequal 
path lengths traversed by the return light 220 from the 
sample volume 213 to the apertures 221a and 222a of the 35 
collection waveguides 221 and 222. The first portion 228 of 
the return light 220 collected by the aperture 221a travels an 
additional path-length of {y l-y2 ? } through the tissue as 
compared to the second portion 230 of the return light 220 
which is collected by the aperture 222a. Hence, the first 40 
portion 228 of the return fight 220 suffers more path-length- 
dependent attenuation as compared to the second portion 
230 of the return light 220. Thus, the signals 224a and 225a 
representing the intensity l cl (X) and l c2 (X) at different wave- 
lengths will exhibit differing levels of modulation caused by 45 
the attenuation of the sample 214. 

In addition, the signals 224a and 225a will also exhibit 
wavelength-dependent modulations caused by the instru- 
mental effects. The wavelength-dependent modulations due 
to the instrumental effects can be determined and, in turn, 50 
compensated for by conducting a calibration of the LIFAS 
system 210. System calibration can be performed using light 
from a standard lamp (Quartz Halogen Lamp, Model No. 
63358, Oriel Instruments, Stratford, Conn.) having a prede- 
termined continuous spectrum to measure the wavelength- 5 s 
dependent instrumental effects of the LIFAS system 210. 
Once these instrumental effects are known, the processor 
223 can be adapted to correct the measured intensities l cl (X) 
and l c2 (X) for modulations caused by the wavelength- 
dependent instrumental effects. The corrected intensities 60 
l ci (X) c and l c2 (k) c representing the intensity of the first and 
second portions 228 and 230 of the return light 220 at 
different wavelengths can then used to determine the attenu- 
ation coefficient a(X) of the sample 214 as described below. 

The regions (T) and (2), of FIG. 2(a) may be at an 65 
effective range, symbolized as "R," from which the collec- 
tion waveguides 221 and 222 collect the majority of the first 



and second portions of return light 228 and 230, respec- 
tively. The effective range "R" will vary with the attenuation 
of the sample. However, the regions ® and (5), are geo- 
metrically symmetric with respect to the aperture 215a of 
the excitation waveguide 215 and, hence, have identical 
intensity of the intrinsic fluorescence 229 represented by 
IjiX). Furthermore, since the lateral and axial distances xl\ 
x2' and y l\ y2* are predetermined by the configuration of the 
probe (in the preferred embodiment xl'=x2', y V>y2\ y2'=0). 
the wavelength-dependent intensities l c] (X) c and l c2 (XY can 
be described by the following equations: 



I fj (Af=MX)^ A,, ' ^+,l, 



11) 

(2) 



The attenuation coefficient a{X) can be calculated inde- 
pendent of the value of t4 R" from (1) and (2) as follows: 



a(^)={l/{y2'-yl'}}ln{I cl (X)^, 2 (?i) c -} 



Where the difference {y2-yl'} is a known constant, V 
is the natural exponential function and "In" is the natural 
logarithm. It will be appreciated by those of ordinary skill in 
the art that the natural exponential and logarithm can be 
replaced by the common exponential and logarithm to the 
base 10, respectively. Similarly, for the probe configuration 
in FIG. 2(c), the attenuation coefficient can be calculated as 
follows: 



(4) 



Once the attenuation coefficient a(X) is determined, the 
intrinsic fluorescence, I/U can be restored from either of 
the signals l cl (X) c or l c2 (X) c (preferably l (2 (XY where 
y2'«yD by assuming an average effective range "R" from 
which most of the intrinsic fluorescence is collected. For 
biological tissue, the constant "R" is approximately about 
0.2 mm at 308 nm excitation radiation. Therefore, the 
intrinsic fluorescence \j(X) can be obtained by substituting 
the measured l c2 (X) c mt0 equation (2) and solving for IjiX) 
using the constant R and the known value y2': 



(5) 



Whether the measured attenuation accounts for absorp- 
tion and/or scattering is mainly determined by the wave- 
length band of interest and the nature of the sample. The 
optical properties of biological tissue and the effects of tissue 
on LIF are greatly different for the wavelengths below and 
above approximately 600 nm. Below about 600 nm. the 
optical attenuation of biological tissue is primarily due to 
absorption and, hence, the attenuation coefficient a(X) will 
represent absorptivity a(X). Absorptivity is a property of a 
substance, while absorbance is a property of a particular 
sample of a substance. Therefore, absorbance will vary with 
the concentration of the substance (e.g. hemoglobin* and 
geometry of the tip of the probe. Thus, for samples such as 
biological tissue where the attenuation of the sample is 
primarily due to absorption, the absorbance A(A) and the 
percent transmittance %T(?i) of the sample can be calculated 
as follows: 



I6l 
(7) 



In an alternative embodiment of the LIFAS system 310, 
shown in FIG. 3(a), the source 311 emits laser radiation 312 
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at a wavelength and intensity capable of inducing fluores- 
cence of the sample 314. The laser radiation 312 is reflected 
by a dichroic mirror 326 into an excitation-collection 
waveguide 321, which transmits the laser radiation 312 to 
the sample volume 313 where the laser radiation 312 excites 5 
local fluorophores to emit intrinsic fluorescence 329. 
Furthermore, the intensity of the laser 312 should be suffi- 
cient to create a sample volume 313 that overlaps with the 
numerical aperture 319 of the collection-only waveguide 
322. The dimensions of the sample volume 313 will also 10 
depend on the numerical aperture of the excitation- 
collection waveguide 321, on the wavelength and intensity 
of the excitation radiation 312 and on the optical properties 
of the sample 314. The intrinsic fluorescence 329 is 
modulated, for example, by the absorption and scattering of 15 
the local chromophores (not shown) and scatterers (not 
shown) of the sample 314. 

The excitation-collection waveguide 321 collects a first 
portion 328 of the return light 320 directly from the sample 
volume 313. The first portion 328 of the return light 320 is 20 
transmitted through the dichroic mirror 326 to the first 
sensor 316. The first sensor 316 generates a first signal 324a 
representing the intensity I, r (k) of the first portion 328 of the 
return light 320 at a plurality of wavelengths within a 
predetermined wavelength band. 25 

As shown in FIG. 3(a), the collection-only waveguide 322 
may be positioned such that the aperture 322a is laterally 
displaced from the excitation-collection waveguide 321 by 
the distance x3 and axially displaced from the aperture 321a 
of the excitation-collection waveguide 321 by the range y3. 30 
Preferably, the lateral distance x3 is zero while the axial 
distance y3 is non-zero. Also, the apertures 321a and 322a 
of the waveguides 321 and 322 are preferably positioned in 
close proximity to or in contact with the sample 314 during 
use of the LIFAS system 310, as shown in FIG. 3(a). 35 
Meanwhile, the numerical aperture 319 of the collection- 
only waveguide 322 is selected to include at least a portion 
of the sample volume 313. 

The collection-only waveguide 322 collects a second 
portion 330 of the return light 320 which is transmitted to a 40 
second sensor 318. The second sensor 318 generates a 
second signal 325a representing the intensity l co Ck) of the 
second portion 330 of the return light 320 at a plurality of 
wavelengths preferably within the same wavelength bands 
as the signal 324a generated by the first sensor 316. 45 
Furthermore, a longpass filter is placed in front of the 
aperture of each of the sensors 316 and 318 to selectively 
block backscattered excitation radiation. 

In the embodiment shown in FIG. 3(a), the first and 
second portions 328 and 330 of the return light 320 expe- 50 
rience different attenuation effects due to the unequal path- 
lengths traversed by the return light 320 from the sample 
volume 313 to the apertures 321a and 322a of the 
waveguides 321 and 322, respectively. The second portion 
330 of the return light 220 collected by the aperture 322a 55 
travels an additional path-length through the tissue as com- 
pared to the first portion 328 of the return light 320 which 
is collected by the aperture 321a directly from the sample 
volume 313. Hence, the first portion 328 of the return light 
320 suffers less path-length-dependent attenuation as com- 60 
pared to the second portion 330 of the return light 320. Thus, 
the signals 324a and 325a representing intensity I ir (A) and 
l co (k) at different wavelengths will exhibit different levels of 
modulation caused by the sample 314. 

The signals 324a and 325a will also exhibit wavelength- 65 
dependent modulations caused by the instrumental effects. 
The wavelength-dependent modulations due to instrumental 



effects can be determined and, in turn, compensated for by 
conducting a calibration of the LIFAS system 310. System 
calibration is performed using light from a standard lamp 
(Quartz Halogen Lamp, Model No. 63358, Oriel 
Instruments, Stratford, Conn.) having a predetermined con- 
tinuous spectrum to measure the wavelength-dependent 
instrumental effects of the LIFAS system 310. Once these 
instrumental effects are known, the processor 323 can be 
adapted to correct the measured intensities I vr (/U and l co (X) 
for modulations caused by the wavelength-dependent instru- 
mental effects. The corrected intensities l xc {X) c and I co {kY 
representing the intensity of the first and second portions 
328 and 330 of the return light 320 at different wavelengths 
can then be used by the processor 323 to determine the 
wavelength-dependent attenuation a(X) of the sample 314. 
As discussed in detail, below, once the attenuation is known, 
either of the signals 324a or 325a, preferably 324a, can be 
corrected for the effects of attenuation to restore the intrinsic 
LIF 329 of the fluorophores in the sample volume 313. 

In the preferred embodiment, the source 311 is an Xe — CI 
excimer laser emitting ultraviolet excitation radiation at a 
wavelength of 308 nm. The waveguides 321 and 322 are 
advantageously comprised of optical fibers or optical fiber 
bundles, which can be integrated into a probe for ease of use 
and durability. A 1 .4 mm diameter optical bundle and a 0.4 
mm optical fiber are used as the excitation-collection and the 
collection-only waveguides, respectively. Preferably, the 
collection-only waveguide is a plurality of 0.4 mm optical 
fibers disposed about the periphery of the excitation- 
collection waveguide. The optical bundle and fibers are 
made of fused silica which is transparent to the 308 nm 
ultraviolet radiation. Furthermore, the probe may be config- 
ured within, for example, the shaft of a hypodermic needle, 
so that the in-vivo physiological or pathological properties 
of biological tissue can be determined. Each of the sensors 
316 and 318 is a spectrograph (Model FF250, ARIES Inc., 
Concorde, Mass.) associated with a 1024 element intensified 
photo diode array (PDA) to image the resolved spectrum. 
Each PDA is connected to an optical multichannel analyzer 
(OMA III, EG&G Princeton Applied Research Corporation, 
Princeton, N.J.) which measures the intensity of the light 
spectrum imaged by the PDA and produce the signals 324a 
and 325a. The processor 323 is a personal computer that is 
networked to each of the sensors 316 and 318 via the signal 
paths 324 and 325 to receive the signals 324a and 325a, 
respectively. 

The method of determining the attenuation coefficient 
a(X) and the intrinsic fluorescence of a sample using 
the embodiment of FIG. 3(a) in accordance with the current 
invention is as follows. In the embodiment shown in FIG. 
3(a), **D" represents the effective range, from which the 
collection waveguide 321 collects a majority of the portion 
328 of the return light 320. The portion 328 of the return 
light 320 is collected directly from the sample volume 313 
by the excitation-collection waveguide 321. Therefore, the 
wavelength-dependent intensity \ KC (Xf of the portion 328 of 
the return light 320 collected by the excitation-collection 
waveguide 321 will be substantially similar to the 
wavelength-dependent intensity of the return light 320, 
represented as IjK). However, the return light 320 with the 
intensity IJX) will suffer additional wavelength-dependent 
attenuation as it travels the extra path-length to reach the 
aperture 322a and is collected as the second portion 330 of 
the return light 320 with the intensity l co (M . The lateral and 
axial distances x3 and y3 will be predetermined by the 
configuration of the waveguides 321 and 322. In the pre- 
ferred embodiment, the lateral distance x3=0 and the axial 
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distance y3>0. Thus, the wavelength-dependent intensity 
can be approximated by the following equations: 



*i tc afe- alXus 

ff(A)=|l/vi}ln|/,,(A)7UAr) 



(8) 



(9) 



(10) 



10 



Similarly, for the alternative probe configuration shown in 
FIG. 3(b), where the lateral distance x3'>0 and the axial 
distance y3'=0, the attenuation coefficient can be approxi- 
mated as follows; 



a(A)={ lA3')MI, r af 



%Ta)^{i co af/ijxf} ioo 



(14) 



15 



(in 



Once the attenuation coefficient a(X) is determined, the 
intrinsic fluorescence, can be restored from either of 
the signals l xc CK) c or l co {XY\ preferably l xr (X) c , by assuming -® 
an average effective range "D" from which most of the 
intrinsic fluorescence is collected. For biological tissue, the 
constant "D" is approximately about 0.2 mm at 308 nm 
excitation radiation. By substituting the measured l xr (X) c 
into equation (7) and solving for I/X): 25 



(12) 



As discussed in connection with the previous 
embodiment, the measured attenuation may account for 30 
absorption and/or scattering depending on the wavelength 
band of interest and the nature of the sample. Below about 
600 nm, the optical attenuation of biological tissue is pri- 
marily due to absorption and, hence, the attenuation coef- 
ficient a(k) will represent absorptivity a(k). Also, the absor- 35 
bance A(k) and the percent transmittance %T(k) of such 
samples can be calculated as follows: 



fJ3) 



40 



In another alternative embodiment shown in FIG. 4, a 
LIFAS system 410 in accordance with the present invention 
has been adapted for biomedical applications, in this par- 
ticular embodiment, the LIFAS system 410 has been adapted 45 
to determine the optical attenuation of biological tissue. The 
laser 411 emits radiation 412 at a wavelength capable of 
exciting the tissue 414 to emit fluorescence. The radiation 
412 is directed through an iris 440 at a dichroic mirror 442 
which reflects the radiation 412 onto a lens 444 which 50 
focuses the radiation 412 onto the proximal tip 445 of an 
optical fiber 446. The adjustable iris 440 is advantageously 
used to reduce the energy of the radiation 412. Alternatively, 
the adjustable iris 440 can be replaced by any suitable 
attenuator. 55 

The probe 448 includes a central optical fiber 446 that is 
used as the excitation-collection fiber and peripheral optical 
fibers 450a-/z are used as the collection-only fibers. The 
distal end of the optical fibers 446 and 4S0a~h are incorpo- 
rated into the optical probe 448. The aperture 448a the 60 
optical fiber probe 448 is placed in proximity to the tissue 
414 so that the aperture 446a of the excitation-collection 
fiber 446 and apertures 451a-/? of the collection-only fibers 
45Qa-h are in contact with the tissue 414. As shown in a 
partial perspective view in FIG. Si a), the optical fiber probe 65 
448 includes a central optical fiber 446 and a plurality of 
optical fibers 450a-/i disposed about the periphery of the 
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central optical fiber 446. The apertures 4Sla-h of the distal 
ends of the collection-only fibers 450a-/? are axially dis- 
placed by a small distance y3 with respect to the aperture 
446a of the excitation-collection fiber 446. The return light 
420 collected by apertures 451a-/? of the collection-only 
fibers 450a-/? pool into the aperture 418a of the sensor 418. 
It will be appreciated by those in the art that the optical fibers 
446 and 450a-/? can advantageously be replaced by optical 
fiber bundles in order to obtain greater flexibility and 
durability. For example, as shown in FIG. 5(b), the central 
optical fiber 446 can be replaced with an optical fiber bundle 
454. 

The excitation radiation 412 is transmitted through the 
optical fiber 446 to the tissue 414 to induce intrinsic fluo- 
rescence of the fluorophores of the tissue 414. The intrinsic 
fluorescence is modulated, for example, by the chro- 
mophores and/or scatterers of the tissue 414. A first portion 
428 of the return light 420 is collected by the optical fiber 
446 from the tissue volume that is directly irradiated by the 
excitation radiation 412 and transmitted through the optical 
fiber 446 to a first lens 444 where the return light is directed 
through the dichroic mirror 442 and then focused by a 
second lens 449 onto the aperture 416a of a first sensor 416. 
Similarly, the collection-only optical fibers 450a-/r collect a 
second portion 430 of the return light 420 and transmit the 
second portion 430 of the return light 420 to the second 
sensor 418. Preferably, longpass filters are placed in front of 
the apertures 416a and 418a of the sensors 416 and 418 to 
selectively block backscattered excitation radiation. A first 
signal 424a and a second signal 425a representing the 
intensity of the first and second portions 428 and 430, 
respectively, of the return light 420 are generated by the first 
and second sensors 416 and 418 and transmitted via signal 
paths 424 and 425 to the processor 423. The processor 423 
uses the first and second signals 424a and 425a generated by 
the detectors 416 and 418 to determine the wavelength- 
dependent attenuation of the sample 414 using equations (8 ). 
(9) and (10) as described in the previous embodiment. 

In a particularly preferred embodiment of the LIFAS 
probe 448 shown in FIG. 8, an illumination source 490 
emitting visible light 492 is mounted on the probe 448 to act 
as a spotlight for the operator. Because room illumination 
can contaminate the spectral measurements of the system 
410 by adding background light, use of the illumination 
source 490 will allow the operator to see and accurate] > 
position the probe 448 under low-light conditions. The 
illumination source 490 is configured to illuminate the 
sample at all times, except when the LIFAS system is 
monitoring return light from the sample 414. 

The biological electrical signal arising in nerves or con- 
tractile tissue, such as muscle, is known as the "action 
potential" and is caused by sudden changes in the ion 
conductivity of the cell membrane. The occurrence of an 
action potential in contractile tissue inmates a contraction. 
For example, in cardiac tissue, the action potential propa- 
gates and spreads in a wave-like manner to induce a local 
myocardial contraction wherever it travels. 

An action potential propagating in the tissue local to the 
aperture 448a of the probe 448 can be detected by the 
electrode 464 which is incorporated into the probe 448 as 
shown, for example, in FIG. 6(a). Voltage alternations 
caused by the occurrence of the action potential are picked 
up by the electrode 464 in reference to the common elec- 
trode 466 and are transmitted to the amplifier 460. The 
common electrode 466 provides the electrical ground for the 
amplifier 460 by maintaining contact with the tissue of 
interest 414 or the whole body. As discussed in connection 
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with FIG. 7, below, the amplified action potential 468 is 
transmitted to the processor 423 to trigger the acquisition 
process of the LIFAS system 410 at a preselected phase of 
the tissue contraction or of the cardiac cycle, whichever is 
applicable. 

In an alternative probe arrangement, a plurality of elec- 
trodes or fibers with a conductive coating can be distributed 
circumferentially about the tip of the optical fiber probe 448. 
For example, as shown in FIG. 6(b\ the probe 448 can be 
equipped with three electrodes 470, 472 and 474 arranged in 
a triangular configuration. The action potentials measured 
by each of the electrodes 470, 472 and 474 are amplified 
through separate channels of the amplifier 460 and trans- 
mitted to the processor 423. As illustrated in FIG. 7, the 
processor 423 processes the received signals to determine 
the direction of propagation of the contraction vector 462. In 
particular, the processor detects the phase lead/lag between 
the action potentials collected by the electrodes 470, 472 and 
474 to determine the orientation of the contraction vector 
462 with respect to the location of the electrodes 470, 472 
and 474. 

For example, the contraction vector 462 shown in FIG. 
1(a) is propagating from the tissue site 476a to the site 478a 
and, hence, as shown in FIG. 7(c), the action potential 414a 
arrives before the action potential 412a which in turn arrives 
before the action potential 470a. Moreover, the phase dif- 
ference or time delay between a pair of action potentials 
indicates how the contraction vector 462 is centered between 
the location of the corresponding pair of electrodes. For 
example, as shown in FIG. l{b\ the contraction vector 462 
propagates from the tissue site 4766 to the site 4786. As a 
result, the corresponding time of arrival of action potentials 
4706, 472Z> and 4746, as shown in FIG. 1(d), will vary. The 
processor 423 processes the signals and indicates the direc- 
tion of the propagation of the contraction to the system 
operator. The direction of propagation can be indicated, for 
example, by a circular array 496 of light emitting diodes 
(LED) 498 mounted cirumferentially on the probe 448, as 
shown in FIG. 8. The processor 423 transmits a signal to the 
LED array 496 so that, for example, only the LED element 
pointing in the direction of contraction propagation would 
glow. 

It will be appreciated that the amplifier 460 can be 
replaced by any device that can measure the electrical 
activity of biological tissue such as a differential amplifier, 
an electrocardiogram (ECG), an electromyogram (EMG), an 
electroencephalogram (EEG). depending on the LIFAS 
application. Furthermore, it will be appreciated that optical 
fibers with a metallic or electrically conductive coating can 
be used in place of one or all of the electrodes 464, 470, 472 
or 474 to measure the action potential of the tissue. It should 
be noted that, conventional ECG is generally not suitable for 
triggering the acquisition of LIFS or LIFAS systems, since 
it does not accurately indicate the instantaneous state of 
myocardial contraction at the sample volume. However, 
customary ECG using limb or chest leads can be used to 
trigger data acquisition of the LIFAS system 410 where the 
sample 414 is non-contractile tissue. 

In the LIFAS system 410 shown in FIG. 4, the light source 
411 is preferably a lamp or a laser that emits ultraviolet, 
visible or infrared radiation. In the preferred embodiment, 
the source 411 is an XeCl excimer laser emitting pulses of 
ultraviolet excitation radiation at 308 nm. Where ultraviolet 
radiation is used, it is advisable that the optical components 
used in the acquisition system be made of synthetic quartz 
(fused silica) to ensure maximal ultraviolet transmission and 
minimal instrumental fluorescence. Alternatively, a nitrogen 
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laser, a helium-cadmium laser, a frequency-multiplied laser, 
a solid-state laser, an arc lamp or a light-emitting diode can 
be used as the light source 411. The energy of the excitation 
light 412 is typically between 0.001-10 m Joules. However, 
it will be appreciated that the selected energy level should be 
low enough to avoid tissue ablation and/or photobleaching 
while still being adequate to produce detectable LIF. 

In the preferred embodiment, the sensors 416 and 418 are 
each comprised of a spectrograph (Model FF250, ARIES 
Inc., Concorde, Mass.) associated with a 1024 element 
intensified photo diode array (PDA) detector. An optional 
low fluorescence, long-pass filter (not shown) with a cutoff 
wavelength above 308 nm, preferably 335 nm (Schott 
WG335). is placed before the entrance slit of each spec- 
trograph to selectively block any backscattered excitation 
radiation from reaching the sensors 416 and 418. The 
entrance slit of the spectrograph preferably has a width of 
100 micrometers. The spectrograph uses a 150 lines per 
millimeter diffraction grating to disperse the incoming return 
light 420 into its spectral components. 

The spectrum formed by the spectrograph is imaged by a 
detector, preferably an intensified linear photodiode array 
(Model 1420, EG&G Princeton Applied Research 
Corporation, Princeton, NJ.) facing the output port of the 
spectrograph. The photodiode array generates a plurality of 
electrical signals representing the intensity of the return light 
420 at wavelengths within predetermined wavelength bands. 
Alternatively, the sensors can be constructed of any suitable 
materials, such as individual light-sensitive diodes with 
appropriate band-pass filters for the analysis of spectral 
bands of the return light or an optical spectrum analyzer 
("OS A") for analysis of a broader spectrum. Selection of the 
return light monitoring device will depend on a variety of 
factors, including cost, accuracy, resolution, and whether the 
user is interested in monitoring a single wavelength, a 
wavelength band or an entire spectrum. 

Although in the preferred embodiments shown in FIGS. 4 
and 5, the optical fiber probe includes a central excitation- 
collection optical fiber or optical fiber bundle with a plural- 
ity of collection-only optical fibers disposed around its 
periphery, it will be appreciated by those of ordinary skill in 
the art that the probe can take many forms. Furthermore, it 
will be appreciated that the central optical fiber can be used 
as the collection-only waveguide while all or some of the 
peripheral fibers can be used for excitation-collection or 
collection-only. The latter arrangement is preferred when 
testing highly-attenuating samples to achieve a better signal- 
to-noise ratio. 

FIG. 9 shows three alternative geometrical configurations 
of the collection-only optical fibers 950a-/i about the 
excitation-collection waveguide 946. In these embodiments, 
it is desirable to couple each of the collection-only 
waveguides 950a-/f to a separate sensor so that the intensity 
of each portion of the return light collected by each fiber can 
be monitored. In FIG. 9(a), the collection-only waveguides 
9$0a~h are arranged about the excitation-collection 
waveguide 946 so that their apertures have a helical con- 
figuration. Because the aperture of each successive 
collection-only waveguide is shifted axially away from the 
aperture 947 of the excitation-collection waveguide 946, the 
return light collected by each of the collection-only 
waveguides 950a-/? will be attenuated in varying degrees. 
Use of a probe arrangement having a plurality of collection 
distances will be useful in measuring the polarization and/or 
the attenuation of especially, a sample having a higher 
sensitivity to attenuation due to absorption. Furthermore, a 
probe incorporating a plurality of collection distances have 



6,124,597 

17 18 

application to a larger variety of samples. For example, the absorbance or percent transmittance can also be utilized, 

where the sample is highly attenuating, the collection-only It has been found that, although the LIFA of ischemic tissue 

waveguides having apertures close to the excitation site can is lower than the LIFA of normal tissue over all wavelengths, 

be used to collect the return light. Where the sample is their difference is most prominent between 350 and 450 nm. 

lightly attenuating, the return light collected by the 5 The LIFAS-derived LIFA spectra of normal and ischemic 

waveguides having apertures that are further from the exci- rabbit kidney are shown in FIG. 12. The LIFA spectra shown 

tation site will be useful in determining the attenuation. in FIG. 12 are acquired using a LIFAS system employing 

In FIG. 9(b), the collection-only waveguides 950a~c are 308 nm excitation radiation produced by an XeCl excimer 
displaced laterally with respect to the excitation-collection laser. This LIFAS system uses a 335 nm longpass filter 
waveguide 946. Such a probe configuration will be useful in 10 (Schott WG335) to cutoff backscattered excitation radiation 
measuring the polarization and/or the attenuation of the form the collected return light. Hence, for this particular 
sample where the sample has a higher sensitivity to attenu- LIFAS system the LIFA values below 350 nm are not 
ation due to scattering. The return light collected by the reliable. For biological tissue, the LIFA values in the wave- 
waveguide 950c will be useful in measuring attenuation of length band about 480 nm have the highest signal-to-noise 
a lightly attenuating sample. On the other hand, the return 15 ratio and hence measurement accuracy. As shown in FIG. 
light measured by the closest waveguide 950a will be useful 12, the LIFA of ischemic tissue is lower than the LIFA of 
in measuring the attenuation of a heavily-attenuating normal tissue over the entire spectrum and is particularly 
sample. low in the region from 350 to 450 nm. 

In FIG. 9(c), the collection-only waveguides 95Qa-c are As demonstrated in FIG. 12, the LIFA, absorbance or 

both axially and laterally displaced with respect to the 20 percent transmittance at a predetermined wavelength or 

excitation-collection waveguide 946. This configuration wavelength bands can be used for the detection of ischemia 

combines the advantages of both of the aforementioned or hypoxia. Furthermore, predictive models, spectral recog- 

configurations shown in FIGS. 9(a) and 9(b). nition techniques and associated classifiers can be applied to 

In an alternative embodiment of the LIFAS system of identify whether a given LIFA, absorbance or percent trans- 

FIG. 4, shown in FIG. 10, a single sensor 1016 and an 25 mittance spectrum has been acquired from normal, ischemic 

associated optical multiplexer 1080 are utilized. The multi- or hypoxic tissue. The classifiers can be initially trained with 

plexer 1080 is used to switch the input of the sensor 1016 LIFA, at a predetermined wavelength or wavelength bands 

between the excitation-collection optical path (1046, 1044, acquired from tissue with a known state of perfusion or 

1042, 1049, 1046') and the collection-only optical path oxygenation. 

(1050) such that the return light from each pathway is 30 Another parameter found to be useful in the classification 

measured sequentially. The processor 1023 sends trigger of normal, ischemic and hypoxic tissue, is the intensity of 

signals to the source 1011, the sensor 1016 and the optical return light from a pair of LIFAS collection pathways (i.e. 

multiplexer 1080 to synchronize their actions. If the source [l cl (k) c , l c2 (k) c ] or [I tr (X) r , l co (h) c ] for the LIFAS embodi- 

1011 is a pulsed laser, two or more radiation pulses are ments described above.) For example, FIG. 13 shows the 

typically required to acquire a single attenuation measure- 35 [I xr (X) c , l co (XY] pair at X=480 nm, symbolized as [1 xr (480f, 

ment using the multiplexed system 1010. A single laser pulse I rt7 (480) r ], measured from hyper-oxygenated (x), normal (o) 

can be used if me laser pulse lasts long enough to sustain the and oxygen deficient (+) rabbit kidney. The [I vr (480) r , 

emission of return light 1020 during the minimum period of l co (4S0) c ] from normal and oxygen deficient tissue tend to 

time required by the optical multiplexer 1080 and associated cluster in two linearly separable regions of the two dimen- 

sensor 1016 to acquire the spectral measurements from both 40 sional I vr (^) r -l c ^(W c space. Thus, a simple linear or non- 

the excitation-collection and the collection only optical linear classifier function can be trained on a set of [l Xf A)c, 

paths. The preferred embodiment of the LIFAS system 1010 l co (k) c ] pairs measured using a LIFAS system from normal, 

is similar to that of 410, except that the system 1010 utilizes ischemic and hypoxic tissue. Other classifiers such as arti- 

an electromagnetic dual port shutter as the optical multi- ficial neural networks (ANN) are also being used. The 

plexer 1080 associated with a single sensor 1016 in a 45 trained classifier function can then be used to classify an 

configuration where the portion 1028 and the portion 1030 unknown [I rr (X) r , l co Ck) c ] pair as normal, ischemic or 

of the return light 1020 are measured sequentially. hypoxic. A "nearest neighbor" (NN) classifier has been 

Diagnosis of Hypoxia and Ischemia found to perform satisfactorily. The NN classifier checks the 

The LIFAS methods and devices of the present invention proximity of an unknown [I vr (480) r , I ro (480) r ] pair to clus- 

can be used advantageously for the in- vivo diagnosis of 50 ters of predetermined [I lr (480) r , I ro (480) c ] pairs measured 

hypoxia and ischemia of biological tissue. Hypoxia is a from known normal, ischemic and hypoxic tissue. Other 

deficiency in the amount of oxygen reaching the tissue, e.g., classifiers such as artificial neural networks (ANN) can also 

due to pneumonia, whereas, ischemia is a localized reduc- be utilized. For myocardial and renal tissue it has been found 

tion in arterial blood perfusion, e.g., due to a narrowing of that it is preferable to use [I rr (X) r , l co Ck c )] measured at 480 

arteries by spasm or disease. Ischemia can also result from 55 nm, to optimize signal-to-noise ratio. However, [I rr (k) c , 

hemorrhage of an arterial wound or during surgical proce- l co Ckf ] at other single or multiple pre-selected wavelengths 

dure that temporarily interrupts the blood flow to a body can also be used. 

region. The present invention describes new criteria and An additional parameter for the classification of normal, 

methods for the discrimination between normal, ischemic ischemic and hypoxic tissue, is the wavelength of the peak 

and hypoxic biological tissue, including, in particular renal 60 transmittance of the tissue, symbolized hereafter as X mclK .r^ 

and myocardial tissue, as follows. especially in the 450-500 nm band. An alternative to X mAr . r 

Urilizing LIFAS methods and devices, it has been found is the wavelength of the peak l co Ck) c in the 450-500 nm 

that ischemic (or hypoxic) tissue exhibits a lower laser band, symbolized hereafter as X mm _ r<7 . Both \ max ^ T m & 

induced fluorescence attenuation ("LIFA") than normal tis- Kzax-co snift towards shorter wavelengths as the hemoglobin 

sue of the same type. Since LIFA is an absolute quantity, its 65 in the tissue becomes deoxygenated. For example, FIG. 14 

magnitude can be directly compared with predetermined shows l co (Xf spectra that are acquired using the LIFAS 

standard values. Similarly, in absorption-dominant media, system shown in FIG. 10 employing 308 nm excitation 
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radiation produced by and XeCl excimer laser. It should be 
noted that the X niax _ co of normal tissue shifts to a shorter 
wavelength as the tissue becomes hypoxic; whereas the 
Kjo.x-co °f normal tissue shifts to a longer wavelength as the 
tissue becomes hyperoxic. Specifically* X ma ^ co varies 
between about 480 and 500 nm as blood or hemoglobin 
oxygenation varies between deoxygenated to oxygenated, 
respectively. A separation border can be identified at about 
489.5 nm to separate Km-co °f iwnnal & ssue (peaks above 
489.5 nm) from hypoxic/ischemic tissue (peaks below 489.5 
nm). A simple classifier can be trained to identify tissue as 
hypoxic if its X max _ co is below 489.5, and vice versa. The 
degree of hypoxia is determined from the magnitude of the 
shift in k max _ co from the normal value, the smaller the shift 
the subtle the hypoxia. 

The presence of renal or myocardial ischemia can be 
detected from the shape of the main lobe of the common LIF 
(e.g. l xc (k) c ) spectrum in the wavelength band 350-450 nm. 
The common LIF can be acquired via an excitation- 
collection waveguide of a LIFAS or a conventional LIFS 
system. For example, FIGS. 16(a) and (b) show typical LIF 
spectra acquired from normal and ischemic rabbit kidneys at 
an excitation wavelength of 308 nm. As shown in FIG. 
16(b\ the shape of the main lobe of the LIF acquired from 
normal tissue is skewed to the right (a positive skewness 
value). However, the shape of the main lobe of the LIF 
acquired from ischemic tissue is almost symmetric (a very 
small skewness value). 

Skewness can be measured by a single parameter which 
indicates the degree of asymmetry of a lobe or curve around 
its mean value. For a symmetrical lobe the skewness is zero, 
however for an asymmetrical lobe the skewness can be 
either positive or negative depending on the shape. The 
skewness "S" is given by: 



1 



(15) 



(16) 



where x is the intensity of the fluorescence at the wavelength 
X, x is the mean intensity of the lobe, X is the wavelength and 
N is the number of wavelength measurements in the lobe. 
Therefore, tissue ischemia can be detected by monitoring the 
skewness of the main lobe of a common ( i.e. modulated) LIF 
acquired using an excitation -collection fiber from the tissue. 
A zero or negative skewness value indicates that the LIF 
spectra is acquired from ischemic tissue, while a positive 
skewness will indicate normally perfused tissue. The bottom 
of the central spectral valley is considered as the baseline for 
the definition of the main lobe as indicated in FIG. 16. 

The intrinsic LIF of normal biological tissue will have a 
spectrum with a single peak resembling that of the main 
tissue fluorophores, collagen and elastin. However, in the 
event of ischemia or hypoxia the intrinsic LIF develops a 
secondary peak or a hump at about 470-490 nm which is 
associated with the peak fluorescence emission of NADH. 
The appearance of this hump in an intrinsic LIF spectrum is 
a direct indication of oxygen deficiency resulting from 
hypoxia or ischemia. Unfortunately, this NADH-related 
hump is not readily detectable in common LIF spectra 
because of an overlapping spurious hump created by an 
intensity dip associated with peak hemoglobin absorption at 
about 410 nm. For example, in common and intrinsic LIF 
spectrum of normal and ischemic rabbit myocardium shown 



in FIGS. 15(a) and (b\ respectively, the intrinsic LIF 
spectrum of normal kidney does not exhibit the second 
hump created by the spectral valley associated with hemo- 
globin absorption. To the contrary, the intrinsic fluorescence 

5 of hypoxic myocardium exhibits a secondary hump indicat- 
ing high NADH concentration. Accordingly, ischemia or 
hypoxia can be detected from a LIFAS-derived intrinsic LIF 
spectrum by assessing its NADH-related secondary hump 
either visually or through simple qualitative analysis. The 

10 magnitude of the NADH hump will indicate the level of 
oxygen deficiency caused by hypoxia or ischemia while its 
absence indicate normal oxygenation. 
Classification using Chemometric Models 
Current ischemia or hypoxia detection techniques use 

15 predetermined ratios of spectral intensities at wavelengths 
associated with NADH, elastin and collagen to discriminate 
normal from abnormal tissue. However, the magnitude of 
the intensity does not fully characterize spectral changes 
associated with ischemia or hypoxia. Instead, it is desirable 

20 to undertake spectral classification based on the shape of 
wavelength bands or the whole spectrum. 

The introduction of the LIFAS methods and devices of the 
present invention enables the application of a new modeling 
approach for the measurement of ischemia and hypoxia. The 

25 approach relies on methods of multivariate linear regression 
(MLR), principle component regression (PCR) and the 
method of partial least squares (PLS). These modeling 
methods are well-established chemometric techniques used 
in the field of analytical chemistry to estimate the concen- 

30 tration of a material in a mixture from absorption spectra of 
the mixture. The algorithms and theoretical basis for PLS 
predictive modeling can be found in Brereton, R. G. Chemo- 
metrics: Applications of Mathematics and Statistics to Labo- 
ratory Systems, New York: Ellis Horwood, 1990. A basic 

35 overview of the PLS regression can be found in Gerald and 
Kowalski, "Partial Least-Squares Regression: A Tutorial" 
Analytical Chimica Acta 185 (1986): 1-17. 

In the preferred embodiment, the MLR, PCR and the PLS 
models are trained with LIFA or absorbance spectra as their 

40 input and the concentration of both oxygenated and deoxy- 
genated hemoglobin as their outputs. The blood content (i.e. 
quantity of blood) in the tissue is determined by combining 
the concentrations of both oxygenated and deoxygenated 
hemoglobin. The MLR. PCR and the PLS models are trained 

45 with LIFA spectra from normal, hypoxic or ischemic rabbit 
heart and kidney. 

Once a model is trained, it can be used to predict the 
concentrations of oxygenated and deoxygenated hemoglo- 
bin if given a LIFA spectrum measured from tissue with 

50 unknown oxygenation. The presence and level of hypoxia 
can be determined by comparing the latter-predicted con- 
centrations values to those known values previously mea- 
sured from normal tissue. Hypoxic tissue has less oxygen- 
ated and more deoxygenated hemoglobin concentrations 

55 than normal tissue. Similarly, the presence and level of 
ischemia can be determined by combining the predicted 
concentrations to determine the total blood content and then 
comparing the total blood content to known values of blood 
content measured in normal tissue. Ischemic tissue will have 

60 a lower blood content than normally perfused tissue. The 
PLS model was found to have the optimal accuracy in 
detecting the presence and level of ischemia and hypoxia. 

Although LIFA or absorbance spectra are the optimal 
input to the MLR, PCR or PLS models for concentration 

65 prediction, it will be appreciated that modulated LIF, or 
percent transmittance spectra can also be used due to the fact 
that, even though common LIF represent fluorescence 



